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Abstract.  Numerical  schemes  for  the  compressible  Navier-Stokes  equations  (CNSE)  are  constructed  on  the 
basis  of  the  kinetic  equation  for  the  Chapman- Enskog  NS  distribution  function  the  macroscopic  variables 
of  which  satisfy  the  CNSE.  It  is  clarified  from  this  approach  that  the  inclusion  of  the  collision  effect  in  the 
numerical  flux  improves  the  accuracy  of  the  scheme.  Then,  a  practically  higher  order  scheme  for  the  CNSE 
is  derived  and  the  existing  first  order  Boltzmann  scheme  for  the  CNSE  [Chou  S.Y.  and  BaganolT  D.,  Journal 
of  Comput.  Phys.  130,  217-230  (1997)]  is  recovered  as  its  simplified  version.  The  numerical  computation 
is  carried  out  for  the  CNSE  derived  from  the  BGK  equation.  Comparisons  are  made  with  the  standard 
solutions  for  the  CNSE,  the  results  of  Chou-Baganoff  scheme,  and  the  BGK  solutions  for  small  Knudsen 
numbers. 


I  INTRODUCTION 

The  Kinetic  schemes  or  Boltzmann  schemes  for  the  fluiddynamic  equation  systems  are  important  byproducts 
of  kinetic  theory  (see  e.g..  Refs.  [1]-  [7]  and  the  references  therein).  In  particular,  the  Boltzmann  scheme  for  the 
compressible  Euler  equations  (CEE)  is  well  known  and  is  studied  in  the  framework  of  Cauchy  problem  of  the 
collisionless  equation.  The  extension  to  the  case  of  the  compressible  Navier-Stokes  equations  (CNSE)  is  also 
tackled.  For  example,  in  Chou-Baganoff  scheme,  the  local  Maxwellian,  which  is  employed  as  the  initial  data 
in  the  computation  for  the  CEE,  is  replaced  by  the  Chapman- Enskog  Navier-Stokes  distribution  function  and 
the  numerical  flux  is  computed  by  using  the  collisionless  equation  as  before.  It  is  demonstrated  in  the  problem 
of  normal  shock  wave  that  Chou-Baganoff  scheme  yields  valid  solutions  for  the  CNSE.  The  satisfactory  viscous 
boundary-layer  profiles  are  obtained  by  the  Gas  Kinetic  BGK  scheme  proposed  by  Prendergast  and  Xu 
This  scheme  is  not  based  on  the  collisionless  equation  and  the  effect  of  molecular  collision  is  taken  into  account 
in  the  numerical  flux. 

Notwithstanding  the  successful  results,  the  theory  of  Boltzmann  schemes  is  not  very  transparent  because  of 
the  indirect  relationship  between  the  fluiddynamic  equation  system  and  the  kinetic  equation  that  the  scheme 
is  based  on.  If  the  scheme  is  constructed  on  the  basis  of  the  kinetic  equation  the  corresponding  macroscopic 
variables  of  which  satisfy  the  fluiddynamic  equation  system,  the  relation  to  the  fluiddynamic  equation  will 
become  obvious  by  construction  and  the  criteria  of  the  approximation  will  become  clear.  When  its  simplified 
equation  is  employed  in  the  construction  of  the  scheme,  the  clear  information  of  its  intrinsic  error  will  help 
avoid  the  meaningless  approximation  in  the  actual  numerical  computation.  On  the  other  hand,  if  we  regard 
the  Boltzmann  schemes  as  the  first  step  of  the  Boltzmann  solvers  to  analyze  the  actual  rarefied  gas  flows  for 
small  Knudsen  numbers,  the  kinetic  equation  should  be  related  to  the  full  Boltzmann  equation. 

In  the  present  paper,  we  will  study  the  Boltzmann  schemes  for  the  CNSE  following  the  above  strategy. 
We  will  construct  the  scheme  on  the  basis  of  the  kinetic  equation  for  the  Chapman-Enskog  Navier-Stokes 
distribution  function  the  macroscopic  variables  of  which  satisfy  the  CNSE.  In  the  course  of  the  construction, 
the  necessity  of  the  inclusion  of  the  collision  effect  for  the  improvement  of  the  accuracy,  which  is  claimed  in 
Ref.  [7],  will  be  clarified.  Then,  a  practically  second  order  Boltzmann  scheme  for  the  CNSE  will  be  obtained 
and  Chou-Baganoff  scheme  will  be  recovered  as  its  simplified  version.  Further,  as  a  byproduct  of  the  present 
study,  a  new  second  order  scheme  (in  time)  for  the  CEE  will  be  obtained  from  the  kinetic  equation  that  has  a 
direct  relationship  with  the  CEE. 
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II  THEORY 


A  Notation 

The  main  notation  employed  in  the  present  paper  is  as  follows:  L  is  the  reference  length  of  the  system 
under  consideration;  po  is  the  reference  density;  Tq  is  the  reference  temperature;  Iq  is  the  mean  free  path  of 
the  gas  molecules  in  the  equilibrium  state  at  rest  with  the  density  po  and  temperature  To;  e  =  ^/wlo/{2L)■, 
Lxi  is  the  space  coordinate  system,  (2-RTo)^/^^i  is  the  molecular  velocity,  where  R  is  the  specific  gas  constant; 
L{2RTo)^^/^t  is  the  time;  po{^RTo)~^^^f{xi,^i,t)  is  the  velocity  distribution  function;  pop,  (2i?To)^/^Ui,  TqT, 
PqP  (po  =  RpoTo),  PoPij,  andpo(2-RTo)^/^Qi  are  the  density,  flow  velocity,  temperature,  pressure,  stress  tensor, 
and  heat  flow  vector  of  the  gas,  respectively. 


B  Chapman-Enskog  expansion 


The  Chapman-Enskog  expansion  of  the  Boltzmann  equation 


dt 


(1) 


for  small  parameter  e  provides  the  materials  for  the  construction  of  our  Boltzmann  scheme.  We  briefly  review 
this  expansion  according  to  Ref.  [8]. 

The  fluiddynamic  solution  of  the  Boltzmann  equation  (1)  is  assumed  to  be  in  the  form: 


/  =  f{th,Dh,€), 


(2) 


where  D  is  the  abbreviation  of  differential  operators  with  respect  to  Xi  and  the  vector  h  represents  the  macro¬ 
scopic  variables,  i.e.,  (ho,  hi,  hi)={p,  pui,  3pT j2  +  pUj),  which  are  deflned  as  the  products  of  /  and  collision 
invariants  {ipo  =  Ij  lAj  =  2, 3),  and  ^”4  =  ^j) 

h=  f  V’M-  (3) 

Jr^ 

The  solution  (2)  depends  on  t  and  Xi  only  through  h  and  its  derivatives  Dh.  It  is  also  assumed  that  the 
conservation  equation  system  for  Eq.  (1)  is  in  the  form: 

^  =  ^{h,Dh,e).  (4) 

The  /  and  #  are  expanded  into  the  power  series  of  e 

oo  oo 

fc=0  k=0 


The  coefficients  fk  (fc  =  0, 1,  2,  ■  •  • )  are  the  solutions  of  the  following  integral  equations: 

0  =  T(/o,/o), 


(6) 


V  (— # 

^  ^dh  ^ 


dfi 

dDh^ 


dx, 


l-\-m=k—l  l-\-m=k 

The  sequence  of  integral  equations  are  solved  from  the  lowest  order  under  the  constrain 

[  V’/odC  =  h,  [  il^fkdi  =  0  (fc  >  1). 

Jk-^  Jr? 

The  inhomogeneous  term  of  integral  equations  (7)  must  satisfy  the  compatibility  condition 

^k-i+  f  = 

Jjl3  OXi 

which  determines 


(7) 

(8) 

(9) 
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C  Distribution  function  for  CNSE 


The  distribution  function  that  yields  the  Navier-Stokes  stress  tensor  and  heat  flow  vector  is  given  by 


fcNSE  =  /o  +  e/i-  (10) 

The  coefflcient  /o  is  the  local  Maxwellian 

/o=  (;^^exp(-<72),  (11) 

where  C  =  and  Ci  =  (6  —  .  The  coefficient  fi  is  the  solution  of  the  linear  integral  equation 

2J(/o,/i)  =  [2{C,Cj  -  (12) 

under  the  constrain  (8).  For  hard-sphere  molecules,  /i  is  given  by 

1  fill  c  Fit 

where  the  functions  A{C)  and  B{C)  are  studied  in  Refs.  [9]  and  [10]  [The  function  v{Q  in  Eqs.  (Al)-(A5)  of 
Ref.  [10]  should  be  2\/2i'{Q].  For  the  BGK  equation 

I =  >»-«' 

/i  is  given  by 

/.  =  -i|2(C.C,  -  (15) 


The  truncated  conservation  equation  system 


dh 

dt 


$0  +  e^i, 


is  the  CNSE  and  its  explicit  form  is 


where 


_a 

dt 


P 

pUi 

PilT  +  uj] 


pUj 

pu^Uj  +  \Pij 

p[|T  -I-  uWuj  -\-  PkjUk  +  Qj 


=  0, 


Pij 


P5,,-7ieT«(^  + 


dui  duj  2  duk 


dxi  dxi 


3  dxk^'’^^' 


Qi 


(16) 


(17) 


(18) 


The  7i  and  a  are  71  =  1.270042427,  72  =  1.922284066,  and  a  =  1/2  for  hard-sphere  molecules  and  71  =  72  = 
0  =  1  for  the  BGK  equation. 


D  Kinetic  equation  for  fcNSE 


Let  us  consider  the  fcNSE  the  macroscopic  variables  of  which  satisfy  the  CNSE  (16).  Making  use  of  Eq.  (7), 
we  find  that  fcNSE  satisfies 


dfcNSE  p  dfcNSE 
dt  *  dxi 


JcNSEifcNSs), 


(19) 
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where 


(20) 


JcNSsifcNSE)  =  +  R, 


•'”  =  £'*  E  +  pi) 

fe=0  l-\-m=k-\-l 

The  functional  form  of  JcNSsifcNSE)  is  given  and  it  depends  on  Xi  and  t  through  the  macroscopic  variables 
for  fcNSE  and  their  spatial  derivatives.  For  general  distribution  function  /  we  can  define  the  operator  Jcnse 
in  the  same  way.  Then,  we  notice  that  fcNSE{xi,^i,t)  is  the  solution  of  the  Cauchy  problem 


f{xi,^i,0)  =  fcNSE{xi,^i.,0).  (23) 

If  the  numerical  scheme  for  the  CNSE  is  based  on  the  above  Cauchy  problem,  its  intrinsic  error  is  zero.  If 
we  employ  its  simplihed  equation,  the  resulting  scheme  has  the  intrinsic  error.  Next  we  show  the  simplified 
equation  together  with  its  intrinsic  error. 

The  solution  of  the  above  Cauchy  problem  is  formally  written  in  the  integral  form  along  its  characteristic 
line; 


=  fcNSEi^i,Xi  —  t^i,0)  + 


{t  -  T)^i,T)dT. 


(24) 


Then,  the  macroscopic  variables  and  their  spatial  derivatives  {k  =  0, 1,  •  •  •  and  V*’  means  identity)  at 
t  =  At  are  given  by 


V^h{xi,At)  =  Yi+Y2 


Yi  =  v'^  f  i}fcNSE{^i,Xi  —  At^i,0)d^ 

^2  =  7^/*  [  [J^  +  .R](Ci,a;i  -  {At  -  T)^i,T)dTd^. 

JR^  Jo 


The  la  is  evaluated  as  follows. 


^2  =  J  Ai)  +  [J^  +  muXi 

=  — +  R]{^i,Xi,0)d^  +  O{At^ 
^  OXi  J  J13 


-  At^i,0))d£  +  O(At3), 

), 


(25) 


(26) 


where  the  trapezoid  rule  and  the  orthogonality  property  tpJ^d^  =  0  and  tpRd^  =  0  are  employed.  The 
above  evaluation  indicates  that  the  influence  of  the  molecular  collision  is  O(At^);  the  intrinsic  error  of  the 
scheme  based  on  the  collisionless  equation  (Chou-Baganoff  scheme)  is  0{At^)  and  thus  it  is  first  order  accurate 
in  time.  If  the  term  is  retained,  the  error  would  be  O(e^Af^).  Furthermore,  if  is  replaced  by  J°,  the 
error  would  be  O(eAt^).  The  third  simplified  equation 


dt  ^^"dxi 


(27) 


yields  the  approximate  solution  of  the  CNSE,  the  accuracy  of  which  is  practically  higher  order  in  time  as  long 
as  g  <  At. 

Finally,  we  remark  on  the  Boltzmann  scheme  for  the  CEE.  In  the  construction  of  the  first  order  scheme 
the  Cauchy  problem  of  collisionless  equation  from  the  local  Maxwellian  is  considered.  Deshpande  and 
Perthame  constructed  the  second  order  Boltzmann  schemes  for  the  CEE  by  the  brilliant  modification  of  the 
initial  data.  On  the  other  hand,  the  local  Maxwellian  /o  the  macroscopic  variables  of  which  satisfy  the  CEE 
satisfies  Eq.  (27)  and  the  intrinsic  error  of  the  scheme  based  on  Eq.  (27)  is  zero,  which  implies  that  the  second 
order  scheme  for  the  CEE  can  also  be  made  by  modifying  the  kinetic  equation  instead  of  the  initial  data. 
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E  NUMERICAL  ART 


Now,  all  the  necessary  materials  are  prepared.  Here,  we  construct  a  Boltzmann  scheme  for  the  CNSE  as  the 
finite-volume  approximation  of  Eq.  (27).  In  order  to  avoid  unessential  complexity,  we  consider  the  spatially 
one-dimensional  case;  the  physical  quantities  are  independent  of  xq,  and  x^-  At  t  =  0,  the  macroscopic  variables 
p{xi,Q),  Ui{xi^Q),  and  T{xi,0)  are  assumed  to  be  piecewise  linear  over  cells  (sj_i/2i •Sj+1/2)  and  the  initial 
data  fcNSE{xi,^i,0)  is  made  accordingly.  Let  the  average  of  h{xi,  t)  over  the  cell  {sj-1/2,  Sj+1/2)  be  denoted 
by  hj{t).  Multiplying  Eq.  (27)  by  ^  and  integrating  the  result  over  the  whole  velocity  space  R^,  over  the  cell 
(■Sj-i/2)  Sj+1/2))  and  over  the  time  interval  (0,  At),  respectively,  we  have 

=  hj{0)  -  '^i^j+1/2  -  Fj-1/2),  (28) 

where  Aa:  =  Sj+1/2  —  Sj-1/2  and  Fj^i/2  are  the  flux  vectors  given  by 

Fj±i/2=  [  I  (29) 

Jo  Jr^ 

We  split  the  flux  vector  Fj^if2  into  two  parts  according  to  the  direction  of  characteristic  line  of  Eq.  (27): 


Fj+1/2  —  +  -^j-Hl/2’ 

E^+1/2  =  /  /  ^  6V’/(s,+1/2 

'  Jo  J^I^O 


T0,^,t)d^dt, 


(30) 


where  Sj+1/2  +  0  and  Sj+1/2  —  0  mean  the  limiting  values;  /  for  Sj+1/2  -f  0  is  computed  in  the  right  cell  and 
that  for  Sj+1/2  —  0  is  done  in  the  left  one.  Recall  that  the  intrinsic  error  of  the  scheme  based  on  Eq.  (27)  is 
O(eAt^).  We  can  employ  the  approximation  of  the  integrand: 


fi^j+1/2  T  0,  ^1,  t)  —  (/o  +  e/i)(sj+i/2  T  0,^,  0)  -|-  t[J°  —  G^~](sj+i/2  T  0,^,  0)  -|-  0{et  + 1^),  (31) 

Then,  we  have 

Pf+1/2  =  ^  -  ^1^1)  ("^'+1/2  T  0,t0)dt  (32) 

The  integration  with  respect  to  ^  can  be  done  in  advance.  In  the  case  of  the  standard  Boltzmann  equation,  it 
is  convenient  to  employ  the  polynomial  approximation  of  the  functions  A{C)  and  B{C)  in  /i.  In  the  case  of  the 
BGK  equation,  these  functions  are  polynomials  of  C  [see  Eq.  (15)].  Although  the  computation  is  tedious  and 
the  resulting  expression  is  lengthy,  however,  the  computer  algebra,  such  as  MATHEMATICA,  is  very  useful 
and  we  are  free  from  the  cumbersome  business. 

For  comparison,  we  present  two  other  formulas  of  the  numerical  flux  without  the  term  J°.  The  first  one  is 

-^f+1/2  =  Ad  j  (/o  -I-  e/i)  (sj+1/2  =F  0,  0)d^,  (33) 

Jii%o 

which  gives  Chou-Baganoff  numerical  flux.  The  second  one  is 

^’,>1/2  =  ^  >  +  1^)  T  0,  C,  0)dC  (34) 

This  is  the  modified  version  of  Chou-Baganoff  numerical  flux  (33),  where  the  slope  in  the  cell  is  taken  into 
account  up  to  the  order  of  At^.  For  e  =  0,  the  formula  (34)  becomes  the  2nd  order  kinetic  flux  vector  splitting 
scheme  for  the  CEE,  which  is  second  order  accurate  in  space  but  is  not  in  time  because  of  the  lack  of  the  term 
J°.  The  formula  (32)  with  e  =  0  gives  the  numerical  flux  for  the  CEE,  which  is  second  order  accurate  in  time 
as  well  as  in  space. 
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FIGURE  1.  Comparison  of  shock  profile  for  M  =  3. 


Ill  NUMERICAL  DEMONSTRATION 

We  consider  the  spatially  one-dimensional  problems  of  normal  shock  wave  and  plane  Couette  flow  and 
carry  out  the  computation  for  the  CNSE  derived  from  the  BGK  equation.  Comparisons  are  made  with  the 
BGK  solutions  as  well  as  the  standard  solution  of  the  CNSE.  In  the  following,  the  physical  quantities  of 
the  gas  is  independent  of  X2  and  X3  and  xi  will  be  denoted  by  x.  In  all  the  computations  of  the  CNSE, 
Aa:(=  Sj+1/2  “  Sj-1/2)  and  At  are,  respectively,  uniform  and  no  slope  limiter  is  employed. 

Figure  1  shows  the  normalized  density  p,  flow  velocity  u,  and  temperature  T  in  the  plane  shock  layer  at  the 
upstream  Mach  number  of  3  (M  =  3).  In  the  figure,  the  reference  length  L  is  equal  to  v^^o/2  (e  =  1)  and 
po  and  To  are  those  at  upstream  condition  {x  =  —00).  The  symbol  indicates  the  result  of  the  present  scheme 
(Aa:  =  0.1,  and  At  =  0.005)  and  the  solid  line  is  the  reference  solution  obtained  by  the  method  employed 
in  Ref.  [11].  The  agreement  is  very  satisfactory  and  the  same  observation  is  made  for  higher  Mach  numbers. 
Needless  to  say,  the  highly  nonequilibrium  flows  is  out  of  the  application  range  of  the  CNSE  and  the  results 
in  the  figure  are  not  physically  correct.  This  figure  shows  the  validity  of  the  present  scheme  as  the  solver  for 
the  CNSE.  The  leading  error  of  the  present  scheme  for  time  step  At  is  O(eAt^)  and  e  corresponds  to  the  local 
Knudsen  numeber,  which  is  0(1)  in  the  problem  of  strong  shock  wave.  Thus,  in  this  problem  the  accuracy  of 
the  present  scheme  becomes  first  order. 

The  accuracy  of  the  present  scheme  (32)  is  practically  higher  order  in  time  as  long  as  e  <  At.  Tables  1-3 
shows  the  result  of  the  Cauchy  problem  of  the  CNSE  with  e  =  0.0001  from  the  initial  condition  p  =  1,  =  0, 

and  T  =  1  +  exp(— x^).  The  time  step  is  given  by  At  =  Ax/8.  For  comparison  the  result  of  Chou-Baganoff 
scheme  (33)  and  that  of  modified  one  (34)  are  tabulated  in  the  tables.  These  tables  indicate  that  the  inclusion 
of  the  collision  effect  improves  the  rate  of  convergence  greatly. 

Finally  we  show  the  results  of  the  Couette  flow  between  two  parallel  plates;  both  plates  are  kept  at  a  uniform 
temperature  Tq;  one  of  the  plates  {x  =  0)  is  at  rest  and  the  other  (x  =  1)  is  moving  at  the  speed  of  1  (the 
dimensional  speed  is  ^/2RTq).  Figure  2  shows  the  results  at  the  Knudsen  number  (=  2e/v^)  of  0.012  (  based 
on  the  distance  between  the  plates,  wall  temperature  To,  and  the  average  density).  In  the  computation  of  the 
CNSE  the  slip  boundary  condition  with  the  slip  coefficients  for  the  BGK  equation  and  the  diffuse  reflection 
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TABLE  1.  The  density  p  at  (x,t)  =  (0.8,2). 


Ax 

Present  [Eq.  (32)] 

Chou-Baganoff  [Eq.(33)] 

Modified  [Eq.  (34)] 

0.2 

0.633401 

0.631299 

0.642383 

0.1 

0.633093 

0.631932 

0.637664 

0.05 

0.633083 

0.632472 

0.635384 

0.025 

0.633083 

0.632769 

0.634237 

TABLE  2.  The  flow  velocity  ui  at  {x,t)  =  (0.8,2). 


Ax 

Present  [Eq.  (32)] 

Chou-Baganoff  [Eq.(33)] 

Modified  [Eq.  (34)] 

0.2 

0.060357 

0.054316 

0.061260 

0.1 

0.059545 

0.056383 

0.060139 

0.05 

0.059458 

0.057865 

0.059770 

0.025 

0.059455 

0.058658 

0.059612 

TABLE  3.  The  temperature  T  at  {x,t)  =  (0.8,2). 

Ax 

Present  [Eq.  (32)] 

Chou-Baganoff  [Eq.(33)] 

Modified  [Eq.  (34)] 

0.2 

1.918041 

1.911346 

1.898904 

0.1 

1.918886 

1.915558 

1.908888 

0.05 

1.918896 

1.917264 

1.913801 

0.025 

1.918887 

1.918081 

1.916315 

is  employed.  The  symbol  of  black  square  indicates  the  result  of  the  present  scheme  for  {Ax,  At)  =  (0.1, 0.01) 
and  that  of  white  circle  does  that  for  {Ax,  At)  =  (0.02,0.002).  They  are  compared  with  the  solution  of  the 
BGK  equation  under  the  diffuse  reflection  boundary  condition  (solid  line).  It  is  seen  that  the  agreement  with 
the  BGK  solution  is  very  satisfactory  (the  Knudsen  layer  correction  to  the  CNSE  solution  is  not  made  in  the 
figure;  this  layer  is  not  taken  into  account  in  the  computation  of  the  CNSE  since  it  does  not  contribute  to 
the  numerical  flux).  Incidentally,  we  remark  that  similar  observations  are  made  in  the  problems  of  weak  shock 
wave  and  nonlinear  heat  transfer  between  parallel  plates  for  small  Knudsen  numbers. 

In  the  present  computation,  the  CNSE  derived  from  the  BGK  equation  is  employed.  The  extension  of  the 
scheme  to  the  CNSE  derived  from  the  standard  Boltzmann  equation  can  be  done  straightforwardly  and  the 
precise  numerical  computation  of  the  standard  Boltzmann  equation  to  establish  the  reference  solutions,  though 
heavy,  can  be  done  at  least  for  spatially  one-dimensional  case  (see  e.g..  Refs.  [13]  and  [14]).  The  extension  and 
the  numerical  demonstration  are  in  preparation. 


IV  CONCLUDING  REMARKS 

We  have  constructed  the  Boltzmann  schemes  for  the  CNSE  that  follow  the  time  evolution  of  Chapman- 
Enskog  NS  distribution  function.  The  agreement  between  the  CNSE  results  and  the  BGK  solutions  seems  to 
support  the  legitimacy  of  this  expansion  up  to  the  Navier-Stokes  level.  As  noted  in  Ref.  [15],  however,  the 
solution  of  the  CNSE  is  obtained  by  the  correction  to  that  of  the  CEE  in  this  expansion.  In  the  abovementioned 
problems,  i.e.,  the  plane  Couette  flow,  weak  shock,  and  heat  transfer  between  parallel  plates,  the  CEE  does  not 
yields  the  first  order  approximation  (weak  solutions  are  not  considered  in  this  expansion).  These  problems  are 
out  of  the  application  range  of  this  expansion  in  the  strict  scense.  The  legitimacy  of  the  use  of  the  CNSE  in 
the  weak  shock  problem  is  given  in  Refs.  [16]  and  [13]  by  the  systematic  expansion  and  it  is  confirmed  that  the 
CNSE  yields  the  solution  correct  up  to  0{M  —  1).  The  CNSE  under  the  slip  boundary  condition  also  yields 
the  solution  correct  up  to  0(e)  in  the  above  two-surface  problems,  which  is  verified  for  the  first  time  by  the 
systematic  analysis  in  the  framework  of  the  boundary-value  problem. 

The  asymptotic  theory  for  small  Knudsen  numbers  has  made  a  great  progress  in  the  last  three  decades  (see 
Ref.  [15]  and  the  references  therein).  The  ghost  effect  is  one  of  the  most  striking  discoveries;  in  Ref.  [17]  the 
discrepancy  of  the  CNSE  in  the  continuum  limit  is  revealed  by  the  systematic  analysis  based  on  the  Boltzmann 
equation  together  and  its  numerical  demonstration.  The  numerical  schemes  for  continuum  gasdynamic  equation 
systems  are  sometimes  employed  as  the  Boltzmann  solvers  for  small  Knudsen  numbers  under  the  seeming 
consistency  with  the  old  theory,  which  is  not  developed  in  the  framework  of  the  boundary-value  problem.  It 
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FIGURE  2.  Comparison  of  flow  velocity  v  and  that  of  temperature  T  at  Kn=0.012. 


is  dangerous  to  employ  the  Boltzmann  schemes  for  this  purpose  without  paying  attention  to  the  recent  theory 
established  by  the  systematic  analyses. 
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